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^ ■ ABSTRACT 



We use the observed sizes and star formation rates of Lyman-break galaxies 



, (LBGs) to understand their properties expected in the hierarchical scenario 

of galaxy formation. The observed size distribution constrains the masses (or 
circular velocities) of the host haloes of LBGs from below, because small 
Qh' haloes can only host galaxies with small sizes. The observed star-formation 



distribution constrains the masses from above, because the number density 



C/3 i of massive haloes in duty cycle is too low to host all the observed galaxies. 

Assuming a flat CDM model with a cosmological constant (JIa ~ 0.7), we find 
I that consistency with observations requires the circular velocities of LEG host 

^ \ haloes to be in the range 100-300 km s^^. The star formation in these haloes 

is quite efficient; a large fraction of all the gas in them can form stars on a 
time scale of about 10-20 per cent of the Hubble time at redshift of 3. The 
predicted comoving correlation length of these objects is ^ 3/i^^Mpc, and the 
predicted velocity dispersion of their stellar contents is typically 70kms^^. 
It is likely that these LBGs are the progenitors of galactic bulges and faint 
ellipticals. 

The same prescription applied to larger haloes in the CDM cosmogony 
predicts the existence of galaxies with star formation rates ~ lOOOMQyr"^ 
at redshift 2; ~ 3. We explore the possibility of identifying these galaxies to be 
the bright sub-millimeter (sub-mm) sources detected by SCUBA. The model 
predicts that the host haloes of these sub-mm sources are massive, with typical 
circular velocity ~ 350 km s^^. The typical star formation timescalc in these 
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systems is about 10 per cent of the Hubble time at redshift 3, and the comoving 
number density of galaxies (in their duty cycle) is 10^''' 10^^/?,'^Mpc^'\ These 
galaxies arc predicted to be strongly correlated, with a comoving correlation 
length of ^ 7/i^^Mpc. They are also predicted to be strongly correlated with 
LEG population at the same redshift, with a cross-correlation length (4 — 
5)/i^^Mpc. The descendants of the bright sub-mm sources should reside in 
clusters of galaxies at the present time, and it is likely that these objects are 
the progenitors of giant ellipticals. We estimate that about 15 per cent of the 
observed sub-mm background is contributed by this population of objects in 
the redshift range 2.5 3.5; a comparable contribution is made by LBGs in the 
same redshift range. The model predicts the existence of a relatively bright 
and red population of galaxies at z ~ 3, which may be observed in the K-band. 

Key words: galaxies: formation - galaxies: structure - galaxies: halo - galax- 
ies: starburst - galaxies: elliptical 



1 INTRODUCTION 

At present there are two main populations of objects which may be crucial to our under- 
standing of the formation of bright galaxies in the universe. The first is the Lyman-break 
population at redshifts z ~ 3, observed by the Lyman break technique (Steidel, Pettini & 
Hamilton 1995). The second is the population detected by some recent sub-mm detectors 
such as SCUBA on JCMT (Holland et al. 1999; Blain et al. 1999). Galaxies in both popula- 
tions are found to show vigorous star formation activity, and so a large fraction of stars we 
observe today may have been formed in these systems 

The Lyman-break technique is quite efficient in identifying star forming galaxies at z ~ 3. 
At the moment, about 1000 galaxies have been selected by this technique and confirmed by 
foUowup spectroscopic observations (Steidel ct al. 1999b). The sample is selected according 
to a flux limit in the optical i?-band {Rab ^ 25.5) which corresponds to UV in the rest- 
frame. Thus, the selection may imply a lower limit on the star formation rates in the selected 
galaxies. The typical star formation rate is ~ 5OM0yr'"^ (all values are quoted for a flat 
universe with Qm.o = 0.3, Qa = 0.7 and Hq — 70kms~^ Mpc~^), assuming a dust correction 
of a factor of about 5 (e.g. Steidel et al. 1999a,b; Pettini et al. 1999a,b; Adelberger & 
Steidel 2000). The comoving number density of the Lyman-break galaxies (hereafter LBGs) 
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is quite high, ~ 2 x 10~^/t^Mpc~^ (Adelberger et ah 1998), comparable to the number density 
of local bright galaxies. They are also strongly clustered in space, with a correlation length 
comparable to that of normal galaxies in the local universe (Adelberger et al. 1998; Giavalisco 
& Dickinson 2000). With a typical half-light radius of about 1.5/i~^kpc (Giavahsco et al. 
1996; Lowenthal et al. 1997), LBGs are quite small compared to the normal galaxies in the 
local universe. 

Since the Lyman-brcak technique selects galaxies according to their rcst-UV luminosity, 
it may miss dust-enshrouded star-forming galaxies, because the UV photons from young 
stars in such galaxies may be absorbed by dust and re-emitted in the far-infrared. Such 
galaxies can be detected in the sub-mm band, and so the recent discovery of strong sub- 
mm sources by SCUBA on JCMT may provide important clues to the star formation at 
high redshifts (see Blain et al. 1999). Unfortunately, current observations are still unable to 
determine the nature of these sub-mm sources. In particular, spectroscopic data on these 
sources are generally difficult to obtain, and so most of them do not have measured redshifts. 
If these strong sub-mm sources have redshifts similar to the LBG population (2; ~ 3), as 
many people believe, the star formation rates in them may be as high as ~ 1000 Mq yr~^. If 
this is the case, the bright sub-mm sources may be an extension of the LBG population to 
higher star formation rates. Observational support for this comes from the fact that galaxies 
with the highest UV star formation rates are indeed detected statistically in the SCUBA 
map (Peacock et al. 2000). 

Theoretically, the nature of the LBG population has been studied by many authors, 
but the results are still inconclusive. There have been suggestions that LBGs might be 
associated with massive dark haloes at z ~ 3 (e.g. Mo & Fukugita 1996; Baugh et al. 
1998). Such association allows one to make specific predictions for the clustering of the 
LBG population, and the model predictions are generally consistent with the observed high 
correlation amplitude (Mo & Fukugita 1996; Baugh et al. 1998; Jing & Suto 1998; Moscardini 
et al. 1998; Haehnelt, Natarajan & Rees 1998; Peacock et al. 1998; Wechsler et al. 1998; 
Coles et al. 1998; Mo, Mao & White 1999; Arnouts et al. 1999; Kauffmann et al. 1999; 
Katz, Hernquist & Weinberg 1999). The star formation properties and structure of the 
LBG population have also been studied under this hypothesis (e.g. Baugh et al. 1999; Mo, 
Mao & White 1999; Shu 2000). However, because of the uncertainties in the observations, 
the theoretical interpretation is not unique. For example, some authors suggest that LBGs 
are starbursts triggered by galaxy interactions (Kolatt et al. 1999; Somerville et al. 2000; 
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Wechsler et al. 2000). In this scenario, the freedom in model construction is larger. For 
example, in one version, all LBGs are assumed to be starbursts triggered by galaxy-galaxy 
collisions, and so the star formation rate is not closely related to the mass of the host halo 
(Kolatt et al. 1999). Alternatively, LBGs may correspond to massive systems which were 
undergoing temporary mergers (either major or minor), as suggested by Somerville et al. 
(2000). 

Theoretical interpretation of the sub-mm sources is even more uncertain. At the moment, 
the observational data are still too sparse to rule out any speculations. It is possible that 
these sources belong to the extension of the LEG population to higher star formation rates, 
as mentioned above, but current observations cannot rule out the possibility that they form 
a completely distinct population. It is therefore important to make model predictions and 
outline key future observations which can be used to distinguish different scenarios. 

In this paper, we use recent observational results to set constraints on theoretical models 
for the LEG population. Our approach is different from those in earlier analyses in that 
it is largely empirical. We show that important constraints can already be obtained from 
current observational data. We use models motivated by the existing observational results 
to make further predictions which can be tested with future observations. In particular, we 
examine the properties of bright sub-mm sources in the same theoretical framework, and 
make predictions for this largely unexplored population. 

The layout of this paper is as follows. We describe our empirical approach in Section 2. 
Predictions for the LEG population and for the sub-mm sources are presented in Sections 3 
and 4, respectively. We summarize and make further discussions of our results in Section 5. 

Since we are dealing with observations at high redshifts, a cosmological model has to 
be assumed in order to transform observables into physical quantities. Throughout this 
paper wc assume a flat model with matter density Om,o = 0.3 and a cosmological constant 
corresponding to fl\ = 0.7. We write the Hubble constant as Hq = 100/ikms~^Mpc~^. 

2 AN EMPIRICAL APPROACH TO THE LYMAN-BREAK POPULATION 
2.1 Observational Inputs 

There are two important observational constraints on the properties of the LEG population. 
The first is their size distribution and the second is their star formation rates. Eased on the 
deep imaging data from the Hubble Space Telescope, the effective radii have been determined 
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Figure 1. The histogram shows tlic observed distribution of the half-light radius of LBGs obtained from the data of Giavalisco 
et al. (1996) and Lowenthal et al. (1997) using a flat cosmogy with Qui,o = 0.3, Qa = 0.7. The solid curve is a log-normal fit 
to the data. 



for 45 LBGs brighter than R = 25.5 (Giavahsco et al. 1996; Lowenthal et al. 1997). In Figure 
1 we plot the observed size distribution, together with a log-normal fit of the data. The 
number of LBGs with determined sizes is still quite small, and so significant revision of the 
observational size distribution is possible in the future. 

The UV luminosity function of the LBG population is straightforward to derive from the 
observed magnitudes (Steidel et al. 1999a,b) and the result is shown in Figure 2. However, 
in order to transform the UV luminosity function into a distribution in star formation rate, 
a correction for dust extinction has to be made. This turns out to be difficult, because the 
correction depends on the details of the dust distribution in individual galaxies. Recently, 
Adelberger & Steidel (2000) carried out a detailed analysis of dust extinction in LBGs, and 
their dust-corrected luminosity function is shown in Figure 2. The average dust correction 
is a factor of 5, compared to the factor of 2 to 3 obtained earlier (Steidel et al. 1998; 
Dickinson 1998). In our analysis, we will adopt the dust correction given by Adelberger & 
Steidel (2000). We use only LBGs with observed Rab magnitude brighter than 25.5. For 
comparison, we also show results for two other dust corrections: a factor of 3 (Dickinson 
1998) and a factor of 11 (Cowie & Barger 1999). For convenience, we will refer these models 
of dust correlation as 

© 1999 RAS, MNRAS 000, 000-000 



6 



0.01 - 



(8 



m,0 A 



dust 
corrected 



0.001 



o 
a 



0.0001 



observed 





Figure 2. The observed and dust-corrected luminosity functions for Lyman-brcak galaxies, adopted from Adclberger & Steidel 
(2000). The solid curves are fits to the data. A cosmology with f2m,o = 0.3 and Qa = 0.7 is used in the calculation. 

• Case (a): a dust correction of a factor of 3; 

• Case (b): the dust correction based on Adelberger & Steidel (2000); 

• Case (c): a dust correction of a factor of 11. 

Once a dust model is adopted, the observed UV luminosity function can easily be converted 
into a distribution of the star formation rate. We use the model as given in Steidel et al. 
(1999b). 

2.2 The Empirical Approach 

Given the size and star formation rate of a galaxy, we use the following procedure to predict 
the other properties of the galaxy. We assume that the star formation in LBGs obeys a law 
similar to the Schmidt relation for local star-forming galaxies: 

EsFR = 2.5 X 10-^° [j^^] ' ' ^oyr^'pc^' > (1) 
where Ssfr is the star formation rate (SFR) per unit area and Egas is the surface density 
of cold gas (Kennicutt 1998). This empirical relation is found to hold reasonably well for a 
variety of star-forming galaxies: from quiescent disks to intensive starbursts. Although the 
validity of this relation for high-redshift star-forming galaxies has not yet been checked by 
observations, there is no obvious reason why star formation at 2; ~ 3 should proceed in 
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a dramatically different way. Assuming a surface-density profile for the gas, we can infer 
the total amount of cold gas in a galaxy from its radius and star formation rate. For an 
exponential surface-density profile (which we adopt in our discussion), the total mass in 
cold gas (Mg) is related to the star formation rate and effective radius it!efr (within 
which half of the UV luminosity is contained) as 



VlO"Me/ l^e^j [Wo) ' 

Notice that this relation is not very sensitive to the density profile assumed, and so is a 

direct result of the empirical star formation relation. If we assume the mass in the cold gas 

to be rUg times the mass of the dark halo, then the halo mass is 

Mh = ^. (3) 

The value of mg is not known. If all gas in a galaxy-sized halo settles into the centre as cold 
gas, then mg = flB/^m,o (where is the density parameter of baryons), which is about 0.1 
if we take the cosmic nucleosynthesis value for a cosmology with Q^.o — 0.3 and h — 0.7. 
At 2; = 3, gas cooling is effective in galaxy-sized haloes, and so we expect mg ~ 0.1 unless 
there is a process which can keep the gas from coohng. In our discussion we will treat mg 
free parameter. 

The observed fiat rotation curves of spiral galaxies suggest that the density profiles of 
dark haloes are approximately isothermal at large radii. If we define the mass of a dark halo 
to be that within a radius rh such that the mean density of the halo is about 200 times the 
critical density (as given by spherical collapse model), then the circular velocity and radius 
of the halo can be inferred from 

V, = [10GH{z)M^]'/^ , rh = V,/[10H{z)] , (4) 

where H{z) is the Bubble's constant at redshift z. For a fiat universe with Qra,o — 0.3, 
H{z)/Ho ~ 4 at 2; = 3. The concentration of a galaxy in its halo is defined to be 

Aeflf = , [O) 

where i?d = i?efr/l-68 is the disk scalclength (Mao, Mo & White 1998). Defined in this way, 
Aeff resembles the spin parameter of disk material in dark haloes (see Mo, Mao & White 
1998). Hereafter we refer Aes as the compactness parameter. 

Based on the above discussion, we can construct a Monte Carlo sample of LBGs by ran- 
domly drawing R^s and from their distribution functions. An extra assumption involved 
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Figure 3. Predicted distribution for the circular velocity Vc of halos hosting LBGs. Panels (a)-(c) correspond to the dust 
correction prescriptions of (a)-(c) defined in section 2.1. For each dust correction model, the solid, dotted and dashed lines are 
for rug = 0.12, 0.06 and 0.03, respectively. 

here is that and R^s are not strongly correlated. This assumption is consistent with 
the observational data for galaxies with i?eff-nieasurements. Our conclusion does not change 
significantly if a positive correlation between and it!eff {ReS oc M^-^, as expected if the 
compactness parameter is not correlated with the halo circular velocity) is included. 



3 PREDICTED PROPERTIES FOR THE LBG POPULATION 

As discussed above, our empirical approach enables us to construct a Monte-Carlo sample 
of LBGs that satisfy the observed distributions of size and star formation rate. Prom this 

Monte-Carlo sample, wc can obtain various other properties of the LBG population and 
analyze their distributions. 



3.1 Properties of the host haloes 

In the model outhned above, LBGs are housed in the central parts of dark matter haloes. 

An important characteristic of the LBG population is the typical mass (or circular velocity) 
of their host haloes. Figure 3 shows the predicted distribution of Vc for the host halos of 
LBGs; three panels (a), (b) and (c) correspond to the three prescriptions of dust correction 
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Figure 4. Predicted distributions of the effective spin parameter Aeff for the host galaxies of LBGs. The symbols have the 

same meaning as in Figure 2. 

as defined in Section 2.1. For each dust correction model, the sohd, dotted and dashed 
hnes denote nig of 0.12, 0.06 and 0.03, respectively. One can see that as rUg increases, the 
underlying host halo circular velocities decrease, as expected from Equations (2-4). Figure 
3 also indicates that for fixed rrig, the required circular velocities of LBG hosts increase as 
the dust correction increases from case (a) to case (c). This trend can be easily understood 
as follows: as the dust correction increases, a given observed UV luminosity corresponds 
to a higher intrinsic total SFR, and to produce the higher SFR, the circular velocities of 
underlying haloes have to increase. More quantatively, for rrig — 0.03, the median values 
of Vc are (180, 205, 250) kms"-^ for the three dust corrections cases (a, b, c), respectively. 
For rrig = 0.06 and 0.12, the corresponding values for Vc are (140, 165, 190) kms^^, and 
(110, 130, 160)kms~^, respectively. In our favored dust-correction model (b), the predicted 
median circular velocity is different from that obtained from early analyses. The reasons for 
this difference will be discussed in more detail in Section 5. 

Figure 4 shows the predicted distribution of the compactness parameter as defined in 
equation (5). For the dust correction case (b), the median values for Aeff are 0.024, 0.028 
and 0.033 for rrig = 0.03, 0.06 and 0.12, respectively. For a given disk mass fraction rrig, the 
predicted LBGs are more compact (i.e., with smaller Aefr) as the dust correction increases, 
because only the more compact objects can produce high (intrinsic) SFR while keeping 
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the size distribution unchanged. For a given dust correction, the LBGs are more compact 
for smaller rrig. This is because the predicted SFR drops as rUg decreases, so in order to 
simultaneously satisfy the observed SFR and size distributions, the objects have to become 
more compact. 

The Aeff-distributions for different dust corrections and mg-choices all have log-normal 
shapes with almost identical dispersion (Jhix^g ~ 0.4. This dispersion is consistent with 
that found in numerical simulations for the spin parameter of dark matter haloes (Warren 
et al. 1992; Cole & Lacey 1996; Lemson & Kauffmann 1998). However, we would like to 
emphasize that the Aefr parameter is related to but not equivalent to the conventional spin 
parameter that describes the angular momentum of dark haloes. Even if the initial specific 
angular momentum of the gas is equal to that of the dark matter, the gas may lose angular 
momentum to the dark halo in the settling process. Further, the self-gravity of the disk 
and the dark matter response to baryonic settling in general reduce the disk size and hence 
decrease the Acg parameter (cf. Equation 5; see also Mo, Mao & White 1998). Also we 
caution that the gas settling process is not well understood, so the final stellar component 
may not be completely supported by rotation if the disk gravity becomes important first; 
this may be the case for many of our systems since they have mg > Xes- In such cases, we 
expect that both the angular momentum and the self-gravity of the gas are important in 
determining the final size. Indeed, if a LEG is not supported by angular momentum, then 
the compactness parameter is merely a measure of the concentration of the star-forming gas 
in a dark halo. 

Figure 5 shows the median values for and Aeff as a function of nig. The dotted, solid 
and dashed lines are for the dust corrections of case (a), (b) and (c), respectively. For all 
the dust corrections, the median value of Vc decreases with increasing rUg while the Aeff- 

dependence on rrig has an opposite sense. Again these trends are a result of the fact that we 
have taken the empirical distributions for the SFRs and sizes of the LEG population (see 
Equations 2-4). 

While in Figure 5 we have plotted the parameter nig over a wide range, this parameter 
is in fact limited by two considerations. First, it must be smaller than or equal to the overall 
baryon mass fraction, /b = ^B/^m,o- For the ACDM cosmogony with Qni,o = 0.3, /b = 0.12, 
if we adopt Ob = 0.019/i~^ as found by Buries & Tytler (1998) from the deuterium abun- 
dance in QSO absorption line systems. Second, as we will show in the next two subsections, 
the duty cycle of LEGs and the abundance of dark haloes require nig ^ 0.06. Hence, the 
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Figure 5. The predicted distributions of the median values of Vc and X^s as a function of nig. The dotted, sohd and dashed 
Unes correspond to the dust corrections of case (a), (b) and (c), respectively. The two vertical dot-dashed lines bracket the 
reasonable range of rrig from 0.06 and 0.12. 

typical value of rrig should be in the range 0.06 to 0.12. These two constraints are shown as 
two vertical dot-dashed lines in Figure 5. The range of rrig implies that the median value of 
the circular velocity Vc for LBGs should be between 130 and 165 kms"^. 

3.2 Star formation timescale and duty cycle 

Since LBGs arc objects with high SFR at 2; = 3, an important question is how long can 
the high star formation activity be sustained. In our model, we can easily estimate the star 
formation timescale t^,, defined as, 

where Mg is the total gas mass in Equation (2) and is the SFR. Note that in our empirical 
approach, the star formation timescale for LBGs is independent of the parameter mg, and 
is only a function of the dust-correction factor. 

The predicted distributions of t-^ in units of the Hubble time ^ 1.5h^^ Gyr at rcdshift 3 
are shown in Figure 6, where the three panels arc again for the three cases of dust corrections, 
respectively. This figure shows that, as the dust correction increases, the predicted star 
formation timescales become shorter on average due to increased intrinsic SFR. In all cases. 
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Figure 6. Histograms of the predicted distributions of the star formation timescales (eq. |6| in units of the Hubble time (tn) 
at redshift 3 for the LBG population. The top, middle and bottom panels correspond to the three cases of dust correction (a), 
(b) and (c) defined in Section 2.1, respectively. 

the star formation timescale is of the order of a few tenths of the Hubble time, typically 
~ 0.2h^^ Gyr. This timescale is somewhat longer than the typical timescales of starbursts 
in the local universe, ~ O.lGyr, but much shorter than the Hubble time. This short star 
formation timescale relative to the Hubble time implies that there must be more haloes 
hosting inactive 'LBGs' than the observed number of LBGs; these haloes are not currently 
forming stars actively so that they are not observable as LBGs. The short duty cycle of LBGs 
has important implications for the abundance of the host haloes, which we will discuss next. 

3.3 Constraints from the halo abundance 

As discussed in the previous subsection, the empirically inferred star formation timescales 
for LBGs are only a few tenths of the Hubble time at redshift 3. If we assume the lifetime of 
a halo is typically a Hubble time, then the short duty cycle means that there must be more 
haloes that once hosted LBGs; most of these hosts are dormant at ^ = 3. Accounting for 
this duty cycle, we can determine the comoving number density of LBG haloes as a function 
of Vc from the dust-corrected luminosity function shown in Figure 2. Figure 7(a) shows the 
required number density of LBG haloes as a function of Vc for the dust correction case (b); 
the solid, dotted and dashed lines are for nig = 0.12,0.06 and 0.03, respectively. The thick 
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Figure 7. The inferred comoving number density of halos hosting LBGs as a function of the circular velocity Vc for the dust 
correction case (b). The solid, dotted and dashed lines correspond to nig = 0.12, 0.06 and 0.03, respectively. The thick solid 
curve is the jHcdictcd halo abundance by the Press-Schecter formalism, (a) only accounts for duty cycle, while (b) ciccounts for 
additional gas reservoir that may exist in LBG host haloes (see Section 3.3). 

solid curve shows the halo abundance predicted by the Press- Schechter formahsm (Press & 
Schechter 1974) in the ACDM model (with Qm,o = 0.3, Qa = 0.7, h = 0.7, and with a shape 
parameter F = 0.2 and a normalization ag = 1 for the power spectrum). 

In the above calculation, all the gas has been assumed to have collapsed into the disk and 
forming stars. However, one can imagine an alternative scenario where initially only some 
gas collapsed into the disk; the gas left in the halo then acts as a gas reservoir to prolong 
the star formation timescale. As argued in §2.1, the total fraction of gas mass that can form 
stars must be smaller than the mass fraction in baryons, /b ~ 0.12. Fig. 7(b) shows the 
predicted number of haloes where we have assumed that all baryons are available to fuel the 
star formation. 

Figure 7 strongly suggests that the mass ratio nig cannot be much smaller than 0.1 if 
one considers the duty cycle of individual LBGs but allows no gas reservoir. Otherwise, the 
required halo number density will exceed the predicted number in the ACDM model. Even 
if additional gas reservoir in a halo is considered, one still requires nig ^ 0.06. 

We have carried out the same calculation for a model with = 1, = 0, /i = 0.5, 
F = 0.2 and ug = 0.6. This model requires nig ^ 0.1 in order to predict sufficient number of 
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haloes with Vc > 200 km s ^. The required rrig is shghtly larger than that allowed by cosmic 
nucleosynthesis . 

3.4 Correlation functions 

Based on the predicted distribution of the circular velocity Vc above, we can calculate the 
correlation function for the LBG population using the method suggested by Mo & White 
(1996; see also Mo, Mao & White 1999, hereafter MMW99). The predicted correlation length 
To as a function of the median value of Vc for the LBG population (assuming varying mg) 
is shown in Figure 8, where we have adopted the dust correction case (b) (the result is very 
insensitive to the assumed dust correction). 

As we have discussed §3.1, the median value of Vc for the LBG population is in the 
range of 130-165 km s^^. From Figure 7, we conclude that the correlation length tq is in the 
range 2.6 to 3.4 h'^Mpc for the ACDM model. This value is smaller than the observed value 
obtained by Adelberger et al. (1998), but is consistent with the more recent determination by 
Giavalisco & Dickinson (2000) who found a correlation length (3.2 ± 0.7)/i~^ Mpc (assuming 
the same flat cosmology as we are using here) for a sample of LBGs brighter than = 25.5. 
The predicted correlation length is also smaller than those given by many earlier studies, 
because of the smaller halo circular velocities we predict for the LBG population. 

Using the model given in Mo & White (1996), one can also predict the correlation 
function of the LBG descendants at the present time, assuming that mergers among LBGs 
are not significant. The predicted correlation length is about 6/i~^Mpc, comparable to that 
of normal galaxies in the local universe. 

3.5 Stellar velocity dispersions 

Our models predict the circular velocities of the haloes hosting LBGs, which are not directly 
observable. A more observationally oriented quantity is the stellar velocity dispersion (cTy) 
weighted by the SFR in the LBG population. As argued in MMW99, the dispersions are in 
general smaller than the underlying halo circular velocities as a result of projection effects 

and the fact that the stars only sample the inner rising part of the halo rotation curve. 
Following MMW99, we estimate ay for each LBG by 




(7) 
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Figure 8. The predicted comoving correlation length ro as a function of the median value of circular velocity Vc- The dust 

correction is taken to be case (b). 




Figure 9. The predicted distribution of stellar velocity dispersions for LBGs. The solid, dotted and dashed lines correspond 
to rug = 0.12, 0.06 and 0.03, respectively. The dust correction is taken to be case (b). 
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where Ssfr is the SFR per unit area, is the total SFR, Vrot('") is the rotation speed at 
radius r, i is the inchnation and i = means a face-on disk. We take sin i to be uniformly 
distributed from to 1. The rotation curve is assumed to be produced by a dark matter 
halo with Navarro, Prenk & White (1997) profile; we take the concentration parameter to 
be 10, the median value found in numerical simulations (e.g. Jing 2000). 

The predicted distribution for is shown in Figure 9 for the dust correction of case 
(b). The sohd, dotted and dashed lines are for = 0.12,0.06 and 0.03, respectively. We 
see that the median vulue of ay decreases with increasing nig because the predicted halo 
circular velocities decrease as rrig increases (see Figure 3). As pointed out in Section 3.1, a 
reasonable range of rrig is between 0.06 and 0.12, which corresponds to a predicted median 
value for Uv of about 70kms~^. 

The stellar velocity dispersions of LBGs may be estimated from the widths of nebular 
emission lines. A preliminary analysis of about a dozen objects gives velocity dispersions in 
the range 60 — 120kms~^, with a typical value of about 80kms~^ (see Pettini 2000). Our 
prediction is compatible with this observational result. Clearly the observational sample so 
far is too small to draw a robust conclusion. We also point out that the correspondence 
between the emission-line widths and the stellar velocity dispersions may not be straight- 
forward (see MMW99 for more discussions) . Obviously, a larger sample is needed to firmly 
test our prediction. 

4 PREDICTIONS FOR BRIGHT SUB-MM SOURCES 

As we have shown in the previous section, haloes hosting LBGs are predicted to have circular 
velocities in the range 100 — 300 km s^^ with a median value of ~ 150 km s^^. As one can see 
from Figure 7, the expected halo number density by the ACDM model in this range is com- 
parable to the required density. However, the ACDM model predicts a substantial number 
of haloes with > 300kms~^. A question naturally arises: what are the observational man- 
ifestations of these even larger haloes? Since the SFR is expected to be a steeply increasing 
function of Vc within our theoretical framework, the SFRs in the central parts of these big 
haloes are expected to be much higher. As mentioned in Section 1, if the bright sub- mm 
sources observed by SCUBA indeed have redshifts ~ 3, their SFRs must be ~ 1000 Mq yr~^, 
much higher than the typical SFR for the UV-selected LBGs. It is therefore very tempting 
to make a connection between the bright sub-mm sources and the massive haloes expected 



© 1999 RAS, MNRAS 000, 000-000 



17 

in the CDM cosmogonies. Below we explore in detail the consequences of this connection. 
Throughout this section, we take the mass fraction of cold gas in dark haloes (mg) to be 
0.1, which is within the plausible range discussed in Section 3.1. 

We identify bright sub-mm sources from Monte Carlo simulations as follows. The halo 
circular velocities are drawn from the distribution predicted by the Press- Schecther formal- 
ism for the ACDM model. Each halo is assigned a compactness parameter Acfr (defined in 
equation 5) drawn from a distribution identical to that for the LBG population. The SFR 
in each halo is derived as described in Section 2. We identify the bright sub-mm sources 
as systems which have SFRs higher than that of the LBGs with the highest SFRs. Since 
the brightest UV-selected LBG has R = 20.75 (~ lOOOMoyr"^) after applying the dust 
correction suggested by Adelberger & Steidel (2000), we choose the bright sub-mm sources 
as galaxies with SFR ^ lOOOMgyr""^. This selection procedure is obviously an oversimpli- 
fied approximation. In reality, whether a source can be detected as a bright sub-mm source 
depends not only on the intrinsic SFR, but also on the dust distribution around the galaxy. 
For example, a galaxy with high SFR may be more likely to be observed as a bright sub-mm 
source when it is more compact. Unfortunately, theoretical models of dust emission from 
star forming galaxies are still very uncertain and it is not yet realistic to adopt a selection 
criterion based on first principles (see, e.g., Granato et al. 2000). Because of this, our the- 
oretical predictions are really only for galaxies with the highest SFRs. The connection to 
the bright sub-mm sources can be made only under the assumption that these sources are 
among the objects with the highest SFRs. With this assumption, many of our predictions 
(such as the circular-velocity distribution and correlation functions) are quite insensitive to 
the details of the selection criterion, because they are generic consequences of the high SFRs. 

4.1 Circular velocities and sizes of the sub-mm sources 

The predicted distribution of halo circular velocities for the bright sub-mm sources is plotted 
in Figure 10. The range of Vc is from 270 to 410kms~^ with a median circular velocity of ~ 
350 km s"^. Compared with the predictions for the LBGs (Section 3.1), the haloes of the sub- 
mm sources are predicted to have much higher circular velocities. The predicted distribution 
of effective (half-light) radius is shown in Figure 11. The median radius, ~ 2.5/i~^kpc, is 
about twice as large as the observed median value for the LBG population (see Figure 1). 
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Figure 10. Predicted distribution of circulcir velocities Vc for the haloes that host bright sub-mm sources. 

4.2 Star formation timescale 

Because the bright sub-mm sources all have high SFRs, it is interesting to examine their 
typical star formation timescale. If we define the star formation time in the same way as 
given by Equation (6), the distribution of for the sub-mm sources looks like that shown 
in Figure 12. The median value of the distribution is about 10 per cent of the Hubble time 
at redshift 3, somewhat smaller than that predicted for the LBGs (Section 3.2). 

Taking into account this duty cycle, the comoving number density of such objects at 
2; ~ 3 is predicted to be about 10~^/i^ Mpc~^ in the ACDM model, which is comparable 
to the observational result of the sub-mm sources with Ss5Qiim ^ 10m Jy assuming their 
redshifts 2; ~ 3 (Dunlop 2000). It must be pointed out that the predicted number density 
is quite sensitive to the SFR threshold adopted in the selection of the sub-mm sources, 
because the massive haloes are already on the tail of the halo mass function at 2; = 3. 
However, the prediction that the haloes of bright sub-mm sources are massive is generic 
because, according to the empirical star formation prescription, the SFR strongly depends 
on halo circular velocity (oc V^-'^) but only weakly on the compactness parameter (oc A~g '^). 
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Figure 11. The predicted size distribution for SCUBA sources. 

4.3 The correlation functions 

Using the model of Mo & White (1996), we can calculate the correlation function for the 
strong sub-mm sources at 2; ~ 3. The inferred correlation length is 6.8 h~^Mpc, about a 
factor of two larger than that predicted for the LBG population. Thus, bright sub-mm 
sources, such as the ones observed by SCUBA, should be very strongly clustered in space. 
The corresponding correlation length for the descendants of these sources at the present 
time, which can be calculated using the method described in Mo & White (1996), is about 
10/i~^Mpc. This correlation length is larger than that for present-day field galaxies, but 
is comparable to that for massive elliptical galaxies. This suggests that the strong sub- 
mm sources may be the progenitors of massive elliptical galaxies observed in clusters of 
galaxies today. With the same model we can also estimate the cross-correlation function 
between the LBG population and the sub-mm sources. We find a correlation length of about 
~ 4.5 h~^Mpc for the model in consideration. 

It should be emphasized here that although the predicted number density of the sub-mm 
sources is sensitive to the selection criterion, the predicted correlation lengths are not. This 
is because while the number density of massive haloes at z ~ 3 depends exponentially on 
the halo circular velocity, the correlation amplitudes have a much weaker dependence on the 
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Figure 12. The predicted distribution of star formation timescales for bright sub-mm sources relative to the Hubble time (tn)- 



halo circular velocity. Thus, the predicted correlation lengths are robust predictions which 
can be used to test the model. 

4.4 The stellar velocity dispersions 

As for the LBG population, we also calculate the SFR-weighted velocity dispersions for the 
sub-mm sources. The expected distribution is displayed in Figure 13. The distribution is 
quite broad with a median of ~ 200kms~^, which is much larger than that predicted for 
the LBG population. This larger stellar velocity is a direct consequence of the fact that the 
host haloes for SCUBA sources are predicted to be much more massive. This prediction is 
also fairly insensitive to the selection criterion. 

4.5 The contribution to the 850//m background 

As we have discussed, bright sub-mm sources have high star formation rate and suffer severe 
dust obscuration. Most UV and optical radiation is absorbed by dust and re-radiated in the 
FIR/sub-mm wavelength. It is therefore important for us to estimate the contribution of 
our model sources to the sub-mm background. According to Scoville & Young (1983) and 
Thronson & Telesco (1986), the FIR luminosity is related to the SFR by 
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Figure 13. The predicted distribution of the stellar velocity dispersions for bright sub-mm sources. 
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Using our Monte Carlo sample, we can estimate the total contribution of massive galaxies 
to the 850/im background by integrating the intensities from redshift 2.5 to 3.5 by 



where L^{Vc, Xes) is the monochromatic luminosity at frequency u of an individual sub-mm 
source, n{z) is the comoving number density of sub-mm sources as a function of redshift 
z, dVcom is the differential comoving volume from z to z + dz and dL{z) is the luminosity 
distance for redshift z. We assume the comoving number density for sub-mm sources to be 
constant within the redshift range z = 2.5 — 3.5. For given SFR, the total FIR luminosity can 
be estimated using Equation (8). However, to convert it to the monochromatic luminosity 
(and flux), we need detailed information about the source spectral energy distribution which 
depends on the uncertain dust temperature. We estimate this by using Equation (7a) in 
Adelberger & Steidel (2000) which connects the monochromatic luminosity to the dust 
bolometric luminosity: Lboi,dust = J^L^K,s5o{z); we take /Cgso = 10 and Lboi,dust ~ 1-47Lfir. 
(see Adelberger & Steidel 2000 for details). 

Carrying out the procedure outlined above, we find that the inferred contribution of 
the sub-mm sources to the 850/xm background is 1.2 x lO'^/i^mJy/deg^. Similarly, for the 
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LBGs, we find a contribution of 1.9 x Tlie contribution of LBGs is about 

the same as that of the sub-mm sources because their larger number density compensates 
their lower SFRs. Both the contributions of the sub-mm sources and LBGs are about 15 
per cent of the observed background, ^ 4.4 x lO^mJy/deg^ (Fixsen et al. 1998). The total 
SFR density for sub-mm sources and LBGs are 0.096, and O.lSMeyr-^/i^ Mpc~^ These 
estimates are consistent with previous studies (Blain et al. 1999). We caution, however, that 
this prediction for sub-mm sources, unlike all the other predictions in this section, depends 
quite sensitively on how we select sub-mm sources in the Monte Carlo simulations. 

5 SUMMARY AND DISCUSSION 

In this paper we have adopted an empirical approach to investigate the physical properties of 
LBG population. We have shown that the observed sizes and star formation rates of LBGs, 
combined with the Kennicutt star- formation relation (Kennicutt 1998), directly constrain 
the circular velocity of the host haloes and the concentration of the observed galaxies within 
them. For the adopted ACDM model, we find that observations require the circular velocities 
of LBG host haloes to be in the range 100 — 300 km s~^, and with a median of ~ 150 km s~^. 
This empirically determined range of circular velocities is different from those found in 
previous studies, because we find that star formation in LBGs has a relatively short duty 
cycle; a more detailed comparison is given in in Section 5.1. The star formation in our model 
haloes is quite efficient; almost all the gas in them can form stars on a time scale ~ 3 x 10^ yr, 
much shorter than the Hubble time at 2; = 3. This means that there must be more haloes that 
can host LBGs than wc actually observed; this has important consequences for the abundance 
of the descendants of LBGs, a point which we return to in Section 5.2. The predicted 
comoving correlation length of LBGs is ~ 3 h~^Mpc, and the predicted velocity dispersion 
of their stellar contents is typically 70kms~^. Both these predictions are compatiable with 
observations (Giavahsco & Dickinson 2000; Pettini et al. 1998a; Pettini 2000). 

The same prescription applied to larger haloes predicted in the CDM cosmogony predicts 
the existence of galaxies with star formation rates ~ lOOOMoyr"^ at rcdshift z ~ 3. We 
identify these galaxies to be the bright sub-mm sources detected by SCUBA. The model 
predicts that the host haloes of these sub-mm sources are massive, with typical circular 
velocity ~ 350kms~^ and typical stellar velocity dispersion ~ 200kms~^. Both the halo 
circular velocity and velocity dispersion are about a factor of ~ 2 larger than those for the 
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LBGs. The typical star formation timescale in these systems is about 10 per cent of the 
Hubble time at redshift 3, and the comoving number density of galaxies (in their duty cycle) 
is lO-^-lO-^/i^Mpc-s. 

Since the LBGs and the bright sub-mm sources are both highly dust-obscured SFR 
galaxies, most of their energies are absorbed by dust and re-radiated in sub-mm sources, 
we estimate that they make comparable contributions (about 35 per cent in total) to the 
sub-mm background. 

5.1 Comparisons with previous studies 

In our model, the predicted circular velocities of the host haloes are smaller than those 
found in most previous analytical and numerical studies (see Section 1 for more complete 
references). For example, by matching the number density of massive haloes with the number 
density of LBGs, the haloes selected as LBG hosts have circular velocities ^ 200kms~^ 
(MMW99; Baugh et al. 1998). The difference arises from several factors. First, our approach 
is more empirically based since we directly use the observed size and SFR distributions. 
For given rrig, the inferred Vc only depends on two assumptions, i.e., the star formation 
occurs in an exponential disk and the Kennicutt (1998) SFR relation applies at z — 3. The 
first assumption is consistent with the surface photometry of LBGs (Giavalisco et al. 1996; 
Lowenthal et al. 1997). The second assumption is plausible, since it is satisfied by local 
starburst galaxies. In any case, this assumption, in one form or the other, has been adopted 
in all studies. We therefore beheve that our estimate of Vc rehes less on theoretical modeUing 
and should be reasonably robust. Second, the dust-correction we favor (Adelberger & Steidel 
2000) is about a factor of two or three higher than that used by previous studies; this implies 
larger intrinsic SFRs for the LBG population. Third, due to the increased intrinsic SFRs, 
all gas in the host haloes of LBGs is typically consumed within a few tens of the Hubble 
time ai z = 3. This means that the duty cycle of the LBG phase must be considered, an 
effect largely ignored in previous studies. A direct consequence of the duty cycle is that the 
number densities of haloes that can host LBGs must be larger than that of the observed 
LBGs (see §5.2). To satisfy this increased number density of host haloes, their circular 
velocities have to be correspondingly reduced. The reduced Vc leads to a correlation length 
ro ~ 3 h~^Mpc, which is smaller than that given by previous studies (e.g. Adelberger et al. 
1998) but consistent with the recent determination by Giavalisco & Dickinson (2000). 
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5.2 Descendants of LBGs and sub-mm sources 

The observed number density of LBGs down to R — 25.5 is ~ 2 x 10~^/i^Mpc~^ (Adelberger 
et al. 1998). As the star formation time scale is only about fifteen per cent of the Hubble time 
at z = 3 in our model (see the middle panel in Figure 6), the haloes of possible LEG hosts 
must be about a factor of ~ 6 more abundant than the observed LBGs. The descendants 
of their stellar components would still have number density of about ~ lO^^/i"^ Mpc~^, if 
merging is not dominant, larger than the number density of local galaxies. It is therefore 
possible that some of the LBGs have survived mergers to become the progenitors of the 
stellar bulges in disk galaxies we see today. The sizes and stellar velocity dispersions both 
match the values for the spheroidal components of disk galaxies such as the Milky Way. 
Further the (comoving) correlation length predicted for their descendants is also similar to 
the present-day normal galaxies. Since the star formation time scale of LBGs is ~ 0.2h~^ Gyr 
with typical SFR of about 5OM0yr~^, the total stellar mass formed is ~ IO^^Mq, comparable 
to the mass of stellar bulges. 

Our model predicts that the bright sub-mm sources are hosted by more massive haloes 
with Vc ~ 350kms~^, and their typical stellar velocity dispersion is ~ 200 km s"-*^. These 
kinematical properties are very similar to those found in elliptical galaxies, therefore it is 
natural to make a connection between the bright sub-mm sources and the present-day giant 
ellipticals (see also Dunlop 2000). The predicted number density of these objects is about 
10~^/i^ Mpc~^ (see §4.1). Because the SFRs in these objects are very high (~ lO^Moyr"^), 
the star formation proceeds rapidly, and all gas is consumed in only about ten per cent 
of the Hubble time at 2; 3. The required duty cycle implies that the number density of 
the remnants of sub-mm sources, if they survive today, can be as high as ~ lO^'^/i^ Mpc~^. 
This is within a factor of few of the observed number density of giant ellipticals in the local 
universe. The total stellar mass formed, ~ 10^^ M©, is also comparable to the stellar mass 
seen in giant ellipticals. 

An important question is what triggers the massive star formation in the strong sub-mm 
sources? In the local universe, the most extreme starbursts are found in dusty ultraluminous 
infrared galaxies, a local analogue of the sub-mm sources. Nearly all these galaxies occur 
in merging galaxies (see Sanders & Mirabel 1996 for a review). In Figure 14 we show the 
merging probability for haloes as a function of the circular velocity Vc at redshift 3 for the 
ACDM model based on the extended Press-Schechter formalism (Bond et al. 1991; Bower 
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Figure 14. The halo merging probabilities as a function of Vc. 

1991; Laccy & Cole 1993). Interestingly, one sees that the maximum merging rate occurs 
for Vc 350kms~^, nearly identical to the predicted median value of Vc for bright sub- mm 
sources. This suggests that massive starbursts may be triggered by major mergers, as seen 
in the local universe. Similar to the local ultraluminous IRAS galaxies, some fraction of the 
energy in them may be contributed by a central AGN fueled by the accreted gas (Genzel et 
al. 1998). 

To conclude, we have modelled the LBGs and sub-mm sources using an empirical ap- 
proach. This approach makes a number of robust predictions about the correlation, sizes 
and kinematical properties of bright sub-mm sources. One interesting consequence of the 
model is the existence of a population of inactive haloes at z ~ 3. The number density 
of these dormant haloes is a factor of 6-10 higher than the observed density of LBGs and 
SCUBA sources. The massive ones are expected to contain large amount of relatively old 
stars and can already be seen in the i^-band. It would be very interesting to explore further 
the observational signatures of these objects. 
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